Quantized Anomalous Hall Effect in Two-Dimensional Ferromagnets 
- Quantum Hall Effect from Metal - 
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We study the effect of disorder on the anomalous Hall effect (AHE) in two-dimensional ferro- 
magnets. The topological nature of AHE leads to the integer quantum Hall effect from a metal, 
i.e., the quantization of a xy induced by the localization except for the few extended states carrying 
Chern number. Extensive numerical study on a model reveals that Pruisken's two-parameter scaling 
theory holds even when the system has no gap with the overlapping multibands and without the 
uniform magnetic field. Therefore the condition for the quantized AHE is given only by the Hall 
conductivity a xy without the quantum correction, i.e., \o- xy \ > e 2 /(2h). 

PACS numbers: 72.15.Rn,73.43.-f,75.47.-m,75.70.-i 



The origin of the anomalous Hall effect (AHE) has been 
a subject of extensive controversy for a long term. One is 
based on the band picture with the relativistic spin-orbit 
interaction [jj, while the other is due to the impurity 
scatterings 0. Most of the succeeding theories follows 
the idea that the AHE occurs due to the scattering events 
modified by the spin-orbit interaction, i.e., the skew scat- 
tering or the side jump mechanism Recently several 
authors recognized the topological nature of the AHE 
discussed in Refs. In this formalism, the Hall 

conductivity a xy is given by the Berry phase curvature in 
the momentum (k-) space integrated over the occupied 
states Also there appeared some experimental evi- 
dences supporting it 0- Therefore it is very important 
to study the effect of the scatterings by disorder, which 
makes k ill-defined, to see the topological stability of this 
mechanism for AHE. 

This issue is closely related to the integer quantum 
Hall effect (IQHE) Q but there are several essential dif- 
ferences. Usually the topological stability which guaran- 
tees the quantization of some physical quantity, e.g., o xy , 
has been discussed in the context of the adiabatic contin- 
uation Hi . Therefore it appears that the gaps between 
Landau levels in pure system are needed to start with 
even though the disorder potential eventually buries it. 
In the IQHE system without disorder, the periodic poten- 
tial is irrelevant because the carrier concentration is much 
smaller than unity per atom. Although numerical simu- 
lations |lC( use lattice models, the main concern is put on 
the limit of dispcrsionless Landau levels separated by the 
gaps. In the present case, i.e., in ferromagnetic metals, 
there are multiple bands overlapping without the gaps 
in the density of states. The periodicity of the lattice 
remains unchanged, which prohibits the uniform mag- 
netic field and also gives a large energy dispersion. In 
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the language of the effective magnetic field for electrons, 
it reaches a huge value of the order of ~ 10 4 Tesla, i.e., 
the magnetic cyclotron length is of the order of the lat- 
tice constant, but the net flux is zero when averaged over 
the unit cell. Therefore these two cases belong to quite 
different limits although the symmetries of the systems 
are common, i.e., the unitary class without time-reversal 
nor spin-rotational symmetry 

In this paper we report on an extensive numerical 
study on two-dimensional (2D) models of AHE includ- 
ing the disorder potentials. It is found that the topolog- 
ical nature of AHE leads to a dramatic phenomenon, i.e, 
the ferromagnetic metal turns into an integer quantized 
Hall system by introducing disorder. This is due to the 
topological stability of the Chern numbers carried by the 
extended states which are energetically separated by the 
continuum of the localized states inbetween. Namely the 
localized state can not have a finite Chern number, and 
the integer topological number can not change smoothly, 
i.e., it jumps when it changes. These two facts leads 
to the protection of the extended state carrying a Chern 
number against the weak disorder. The finite-size scaling 
analysis is c omp atible with the two-parameter RG theory 
of Pruisken |ll| , which predicts the plateau transition at 
\vxy\ — 0.5e 2 //i. The critical exponents are consistent 
with that of the IQHE. This problem is not an academic 
one; the recent technology can fabricate very fine thin 
films of ferromagnetic metals with large enough single 
domain. When the coherent length of such a thin film is 
longer than the thickness, it can be regarded as a multi- 
channel 2D system. These systems can offer a possible 
laboratory to test our theory. 

The essence of the AHE is that the Berry phase of the 
Bloch electron is induced by the spin-orbit interaction 
in the presence of the magnetization, which is modeled 
by the complex transfer integrals 5j. Each band often 
gains finite Chern number even though the density of 
states has no gap. The minimal model describing this 
situation is that proposed by Haldane 0] and its ex- 
tension. This model is defined on a honeycomb lattice 



FIG. 1: Haldane's model defined on honeycomb lattice [T^|. 
Open and closed circles represent the A and B sublattice sites 
respectively. The thick arrows 01^2,3 an d are the lat- 

tice vectors of honeycomb lattice and those of triangle sub- 
lattice respectively. The dashed lines represent next-nearest- 
neighbor hopping. 



containing two atoms in a unit cell (Fig. 1). Using the 
Fourier transformation ct = (c -,, c 1 -) of the spinor rep- 

resentation (eL- , , cL. D ), its Hamiltonian is written as 
H 



Y.k c \ n k c k' where 



■tr^o 



2t\ cos cj) cos(fc ■ a. 

i=l,2,3 



2t\ sin tj) sin(A; 



(1) 



t° is the unit matrix and r 1 ' 2,3 are the Pauli matrices. 
Here we assume the perfect spin-polarization and use the 
spinless fermions. The complex next-nearest neighbor 
hopping integral t\e 1 ^ is introduced in addition to the real 
one to between the nearest neighbors. o^ r 2 3 are the lat- 
tice vectors of triangle sub-lattice, while a\ c 2 3 are those 
of honeycomb lattice. 

The extended model is given by adding another layer 
with the change t\ — > —t\ to the original single-layer 
model given above. Furthermore we introduce the energy 
difference between the layers by shifting the uniform po- 
tential ±u. Then the extended model has the symmetric 
and gapless density of states in contrast to the original 
one. In Figs. El are shown the density of states and a xy 
for the single-layer (a) and double-layer (b) models, re- 
spectively. In the single-layer case, a xy is quantized to 
be e 2 /h when the Fermi energy lies within the gap, while 
it is not in the double-layer case where the gap collapses. 

Now we introduce the on-site disorder potential to the 
single-layer model, which is randomly distributed in the 
range [—W/2,W/2], and study the localization problem 
in terms of the transfer matrix method Figure Efa) 
shows the dependence of the renormalized localization 
length Xm/(2M) of a quasi-lD tube with 2M-sites cir- 
cumference on the strength of disorder potential W. In 
each figure, the lines for M = 4,8, 16, 32 are plotted. We 




FIG. 2: Hall conductivity a xy and the density of states for 
the pure single-layer model with t\ = 0.2to, m = 0.4io and 
4> = tt/3 (a), and for the pure double-layer model with t 1 = 
0.2t 0) m = 0.4to, u = 0.8t and <f> = tt/3 (b). 



have also calculated up to the system size M = 64 around 
the lower extended state. (These additional data are used 
in the analyses for the localization length and its critical 
exponent.) The length of a tube is typically ~ 10 5 sites 
and the accuracy of data is within a few percent. We 
can see that extended states are isolated in energy. They 
merge with each other at a critical value 6 < W c < 7, and 
disappear, i.e., all states are localized. This behavior is 
the same as that observed in the ordinary IQH system 
on a square lattice with external magnetic field |H} . It is 
noted that the pair annihilation of extended states always 
occurs between those with the opposite Chern numbers 
[TfH |. Actually, two extended states in Fig. originate 
from lower and upper bands which have opposite Chern 
numbers, ±1, respectively. 

Next we analyze the data to obtain a characteris- 
tic length £(E, W), which depends on E and W but 
not on M, by the scaling hypothesis, Xm(E, W)/M = 
f (£(E, W)/M) , where f(x) is a scaling function [lj. As 
for a localized state, £(E, W) is interpreted as its local- 
ization length in the thermodynamic limit. Figure Etc) 
shows £(E, W = 5.0io) around the lower extended state 
at E — E c . From this data, the critical exponent v 
(£ oc \E - E c \- V ) is estimated as v = 2.37 ± 0.05 with 
E c = (-0.69±0.01)i Q . FigureUd) gives the log-log plot 
of the localization length £(E, b.Oto) as a function of the 
energy measured from E c . The fitting is also shown as 
a solid line with the slope — v = —2.37. This value is in 
reasonable a gree ment with that estimated in the ordinary 
IQH system 13, e.g. v = 2.35 ± 0.03 

The system-size dependence of <r xy shown in Fig. EJb) 
represents its scaling property. There are two critical 
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FIG. 3: (a) Localization length Am of a quasi-lD tube 
where M is the number of A(B) sites on the circumference. 
M — 4, 8, 16, 32 are plotted, (b) System-size dependence of 
a xy in the single-layer system with 2L x 2L lattice points. The 
numbers of samples averaged are 81920, 20480, 5120, 1280 for 
L = 4, 8, 16, 32 respectively. The errors are one standard devi- 
ation, (c) Log plot and (d) Log-log plot of Localization length 
£ around the lower extended state at E c — (—0.69 ± 0.01)£o 
The solid line is a fitting result with the slope — v = —2.37. 



points where there is no size dependence and which sep- 
arate the two energy regions with the opposite size depen- 
dences. a xy at these critical points takes the value about 
0.5e 2 /h. This is consistent with the analysis in terms of 
the effective field theory for the ordinary IQH system in 
the weak-localization region and strongly suggests 
that the critical properties of this transition are same as 
those of the plateau transition a xy : <-> 1. The energy 
of these critical points coincide with that where the local- 
ization length diverges in Figs.OJa) and (c). This means 
that the extended states with Chern number ±1 exist 
there, and a xy in the thermodynamic limit (M — > oo) 
stays quantized to be e 2 /h between these two energies 
[broken line in Fig.OJb)]. 

Now we shall consider the double-layer model. The 
IQHE never occurs in this system with the above set of 
parameters unless the effect of disorder is taken into ac- 
count. We consider the nontrivial case in which there are 
scattering events both within and between these layers. 
In this case, there is no gap between the initial and final 
states of the elastic scattering, and it is possible that all 
the states are localized once the disorder is introduced. 
However as shown below, the extended states and the 
Chern number carried by them are stable against the 
weak disorder. Here we define the strength of intra-layer 
scattering as Wo and represent the strength of inter-layer 
scattering by W\ . As seen in the upper panels of Fig. 0] 
there occur two energies, i.e., E = — 1.5trj and E = 0, 




-3 -2 -1 -3 -2 -1 
E/t n 



FIG. 4: Upper panels : the localization length Am of a quasi- 
lD tube where M is the number of A(B) sites on the cir- 
cumference. M = 4,8, 16, 32 are plotted. Lower panels : the 
system-size dependence of a xy in the double-layer system with 
2L x 2L lattice points. The numbers of samples averaged are 
81920, 20480, 5120 for L = 4, 8, 16 respectively. The color of 
each lines is for the same system-size as in Fig. EI a) and (b) . 



at which Am/2M does not show M-dependence. This 
means that the extended states survives there as in Figs. 
3 at least up to W\ = l.Oio- 

We next present the system-size dependence of a xy , 
as shown in the lower panels of Fig. There appears 



two critical points for the transitions a x 







1 



and a xy : —1 <-> 1. From the particle- hole symme- 
try, a xy {—E) = —<j xy (E) is concluded. Therefore, there 
appear three critical points in the whole energy region. 
Because the single-layer model has at most two critical- 
points, the double-layer model could have maximum four. 
One may wonder why there appear only three critical- 
points in the present case. This is because the middle 
one is composed of two extended states which originally 
contribute to different critical points but carry the same 
Chern number, i.e. —1. These extended states merge 
(at least in the present numerical accuracy) but never 
pair-annihilate, because the composite of these extended 
states carry the non-zero Chern number —2. In other 
words, the conservation law of the topological charge pre- 
vents the localization. 

a xy at the lower critical-point takes the value = 
0.5e 2 //i. This value is again consistent with the anal- 
ysis by Pruiskin and coworkers lllj . However, a xy at 
the middle critical-point is zero. This critical behavior 
seems to violate the prediction by the analysis in [Tl) . 
However, recent numerical studies for the ordinary IQH 
system reveals the new type of critical phenomena, i.e. 
the direct transitions a xy : <-» n (n > 1) |10| . which 
were experimentally observed in advance 17]. The crit- 
ical property around the middle point is considered to 
belong to the same class as a xy : <-» 2. Although the 
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sample size is not large enough in the double-layer model, 
the size-dependence of a xy is consistent with the quan- 
tized plateau shown by the broken line in the lower panels 
of Fig. H 

It is not difficult to generalize the non-linear sigma 
model approach for the localization problem to the case of 
multi-component model without time-reversal nor spin- 
rotational symmetry. This "components" means orbitals, 
spins, and channels in the multilayer cases altogether. 
This approach does not assume the finite gap at the start- 
ing. Following the derivation in Ref. JlJ, we obtain the 
Lagrangian: 



L[{Qi}} 



l.l/ 



+Tr In 



E-H 



irjs 



(2) 



where [g\u> is the scattering strength between compo- 
nents / and I', and [//];';» = Si'iSfi. Tr (In) is the trace 
(logarithm) of matrix with functional index (r) and dis- 
crete indices (p,a,l), where p = ± corresponds to the 
advanced and retarded fields respectively, and a runs 
over replicas. Tr O is the abbreviation for J dr TrO(r) 
where 0(r) is a matrix with p and a indices. The non- 
linear sigma model is the effective model for the massless 
Goldstone modes. In order to extract these modes, the 
parametrization Qi = is useful. From the above 

Lagrangian, it is clear that inter-component scatterings 
Iff" 1 ];;' (I l° c k the out-of-phase modes Ti ^ Ty 
(I ^ I'), and therefore the effective model for massless 
in-phase modes, i.e., Ti = T, reduces to the model iden- 
tical to that in Ref. 11]. The coefficients of the stiffness 
and topological terms for these modes coincide with a xx 
and a xy respectively. It is noted that these a xx and a xy 
contain the contributions from all components, i.e., all 
orbitals, spins and channels. Then the scaling of a xy and 
o~ xx remains the same as given in Ref. . This supports 
the finite-size scaling study given above. 

In real systems, the Coulomb interaction can not be ne- 



glected. In IQHE system, the InT dependence of a xx is 
observed llq | and is attributed to the quantum Coulomb 
correction However, the quantized a xy in the ground 
state is well described by the noninteracting electron 
model. The situation is similar here for the quantized 
AHE system. In the thin film of Fe, In T-dependence of 
a xx is observed while not for a xy 20] , which is explained 
by the quantum Coulomb correction combined with the 
skew scattering mechanism pl|. 



Usually AHE is estimated by pn 



xx u xy, 



where 



Ph, Pxx and a xy are measured as quantities in 3D. In 
good metals, p xx is very small at low temperatures, and 
hence \cr xy \ is large although \ph\ is very small. There- 
fore it is possible that the quantized AHE is realized even 
in the conventional metallic ferromagnets such as Fe or 
Ni, when the thin film is fabricated. Actually, when we 
virtually consider thin film of n-layer systems, the 2D 
\a xy \ at Tc/2 is estimated as ~ 0.59ne 2 /h for Fe p^ . 
- 0A7ne 2 /h for Ni and ~ 0.20ne 2 /h for SrRu0 3 
(from the first article in Ref. 0)- Therefore, the con- 
dition \a X y\ > 0.5e 2 /h is not so difficult to achieve in 
the thin films of metallic ferromagnets. Extrapolating 
the InT behavior of a xx experimentally observed [20| . 
the crossover temperature T croS s from weak to strong lo- 
calization is estimated as T cmss = T e , where To a 
reference temperature of the order of 10K, a® the Drude 
conductivity at To, and A is a sample-dependent scaling- 
exponent of the order of unity. Therefore, if the minimal 
aoh/(Ae 2 ) is less than ~ 10, we have the chance to ob- 
serve the quantized anomalous Hall effect in the exper- 
imentally realizable temperature. Considering that a xy 
have to be larger than 0.5e 2 //i, this condition means that 
the ratio ao/a xy should be smaller than ~ 10. This novel 
quantized Hall state would prove the most dramatic con- 
sequence of the topological nature of the AHE. 
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